Deviation from Snell's Law for Beams Transmitted Near the Critical Angle: 

Application to Microcavity Lasers 
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We show that when a narrow beam is incident upon a dielectric interface near the critical angle 
for total internal reflection it will be transmitted into the far-field with an angular deflection from 
the direction predicted by Snell's Law, due to a phenomenon we call "Fresnel Filtering" . This effect 
can be quite large for the parameter range relevant to dielectric microcavity lasers. 
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A promising approach to making high-Q optical micro- 
cavities and micro-lasers is to base them on totally- 
internally reflected modes of dielectric microstructures. 
This approach is currently under intense investigation 
with resonators and lasers based on a range of shapes: 
disks, cylinders, spheres deformed cylinders and 
spheres g-Q] squares || and hexagons Many dif- 
ferent mode geometries have been observed in such res- 
onators, e.g. whispering gallery modes |1J, bow-tie modes 
H,|J, triangle ||f7| and square modes Typically these 
modes correspond to ray trajectories which are incident 
on the boundary of the resonator above or at the critical 
angle in order to achieve adequately high Q-values and 
may correspond to periodic ray "orbits" (POs) which are 
either stable, unstable or marginally stable. The natu- 
ral and simplest method for predicting how such a mode 
will emit or scatter light is simply to apply Snell's law 
to the underlying ray orbit and follow the refracted ray 
into the far-field. For a ray which is incident at the criti- 
cal angle this would imply emission in the direction tan- 
gent to the emission point. However in several recent 
experiments very large deviations from this simple ex- 
pectation were observed We show below that such 
observations may be explained as arising from the an- 
gular spread in the resonant mode around the PO, and 
the very rapidly varying transmission probability as a 
function of angle near the critical angle. This "filters" 
out the angular components which are totally internally 
reflected (TIR) and preferentially transmits those which 
are far from TIR, leading to a net angular rotation of the 
outgoing radiation from the tangent direction. We call 
this effect Fresnel Filtering (FF). 

The basic effect occurs for a bounded beam of arbi- 
trary cross-section incident from a semi-infinite medium 
of index n into vacuum, although it will be quantitatively 
altered in a resonator due to the curvature and/or finite 
length of the boundary. We thus begin with the infinite 
planar example, which we can solve analytically, before 
presenting numerical results for quadrupolar asymmetric 
resonant cavities (ARCs) There is a large litera- 

ture on reflection of a beam from a dielectric interface 
near or above the critical angle, as the reflected beam 
exhibits the Goos-Hanchen lateral shift as well as other 



"non-specular" phenomena pT| . However only a few of 
these works address the transmission |l^|l3|] of the beam 
and these tend to focus on the evanescent effects in the 
near field; none appear to have identified the Fresnel Fil- 
tering effect and its relevance to dielectric micro-cavity 
resonators. 

For simplicity, we consider a two-dimensional planar 
interface which separates two semi-infinite regions with a 
relative dielectric constant n. Consider a beam Ei a inci- 
dent from the denser region with a central incidence angle 
9i. We will take the beam to be gaussian with a mini- 
mum beam waist w (which we will use to scale all lengths 
henceforth) at a distance z from the interface. The basic 
effect is independent of the nature of the input beam as 
long as it is focused and has significant angular spread. 
The corresponding Snell emission angle 9 e (which is in 
general complex) is given by rising = sin6> e . Si : (xi, zi) 
and S e : (x e , z e ) refer to coordinates tied to the incident 
and refracted beams respectively (see Fig. [j] inset). We 
will consider linearly polarized beams, the correspond- 
ing beam fields E a ,a — TM,TE, will then denote the 
electric (TM) or the magnetic (TE) fields normal to the 
plane of incidence. 

Using the angular spectrum representation [j[4| , the 
incident beam in Si coordinates will consist of a super- 
position of plane waves of the same frequency ui with 
a gaussian distribution of transverse wavevectors nk Q s, 
where s — sin A9i , k Q = to / c a is the wavevector in vac- 
uum and A9i is the deviation angle of the plane wave 
component from 9f. 
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where c = y/i — s 2 and the dimcnsionless width parame- 
ter A = nk w. 

The beam on the z > side of the interface in polar 
coordinates (p, <fi) attached to the interface (after refrac- 
tion) is then given by the integral: 
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exp 



i—pcos (4> — 9 e — A9 e ) 
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Here A9 e is obtained from nsm(8i + A9i) = sm(9 e + A8 e ) 
and Q(s) is given by: 



Q(s) = exp 



— ) s 2 + iAVT- -s-:., 



(3) 



Evaluating this integral in the asymptotic farfield (p — > 
oo ) using the saddle point method we obtain our "gaus- 
sian model" (GM) for FF field: 
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where the transmission functions, evaluated at the rele- 
vant saddle point 
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are given by 



T a {s (4>)] = 



2n\J n 2 — sin 2 cf> 
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sin 2 (j) + n 2 \/l — sin 2 i 
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Here, p = 1 for a = TE and /i = n for a = TAf. 
The relevant saddle point arises from setting to zero the 
derivative of the cosine in the exponent of Eq. (^) ; this 
saddle point value selects the angular component which 
refracts into the observation direction (j> by Snell's law. 
However the amplitude factor obtained by gaussian in- 
tegration around the saddle point shifts the maximum 
of the outgoing beam away from the Snell direction. As 
noted, the effect occurs for narrow beams with an ar- 
bitrary (non-gaussian) wavevector distribution Bis); in 
such a case the factor Q{s ) in Eq. (0) is replaced by 
B(s ) (see e.g. ref. @). 

Eq. (^) gives the angular beam profile in the far- 
field, which is non-zero for any incident angle Oi, even 
9i > 9 C = sin _1 (l/n). The key point is that the angu- 
lar maximum of this outgoing beam, ^maxt *s in gen- 
eral not at the angle 9 e predicted by applying Snell's 
law to the central incident beam direction 9i . Instead, 
due to the unequal transmission of the different angu- 
lar components, the beam direction is shifted by an an- 
gle A9ff corresponding to less refraction than expected 
from Snell's law. This angular deflection can be quite 
large for incidence near 9 C in typical microcavity res- 
onators; in Fig. [l| the dashed line is the result of Eq. (^) 
for critical incidence, for which the Snell angle is <f> = 90°, 



but 4> max = 62° giving A9 C FF = 28°. The farfield peak- 
shift A9ff depends on the beam width A and on n; 
analysis of the stationary phase solution gives the result 
that at 9j = 9 r 



A9 c t 



'ff ~ (2/ tan$ c ) ' A ' (7) 

which predicts A9 C FF w 30° for the parameters of Fig. [|. 

Two technical points are in order here: First, while 
deforming the contour of integration in Eq. (|^) to the 
steepest descent path, depending on how 9 i} 9 C and <j> 
are situated, one might intercept branch cuts due to the 
branchpoints s = ±1 and s = ±sin(#j — 9 C ). Such 
branchpoint contributions are well-studied for the re- 
flected beam shifts [ pi|Jl(|| , but the contribution is sub- 
dominant with respect to the first order asymptotic term 
derived in Eq. (|j); we neglect such terms here. Second, 
there is another saddle point s = cos(9i) which cor- 
responds to angular components with grazing incidence 
to the interface. Because the Fresnel transmission fac- 
tor vanishes for such components s~o only contributes to 
the integral at order 0(p~ 3 / 2 ). This contribution only 
becomes important very near cf> — n/2, where the dom- 
inant contribution from so also vanishes, and again we 
neglect such terms in the current work. 




FIG. 1. Angular farfield intensity distributions 
I ((f)) — \E e ((j>)\ 2 for: (dotted line): critical incidence on a 
planar interface with n = 1.56, A = 8.82, zo = 5 using the 
gaussian model (Eq. (Q)). (solid line): Exact quasi-normal 
mode with diamond geometry at nk r m 90, for a quadrupo- 
lar ARC with e = 0.1, n = 1.56. (dot-dashed line): Chiral 
version of diamond resonance (see text) which eliminates in- 
terference effects. Inset: Coordinates and variables for the 
GM calculation. 

Clearly the same Fresnel Filtering effect will occur in 
emission from dielectric resonators with a magnitude sim- 
ilar to the planar case when the typical radius of curva- 
ture is much larger than w. As an example, we inves- 
tigate the effect of FF on the farfield emission pattern 
of quadrupolar ARCs , dielectric cylinders with 

cross-section given by r((f> w ) = r e (l + £COs2^> w ). We 
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study the exact numerically generated quasibound TM 
modes of a resonator with 10% (e = 0.1) deformation 
for different values of the refractive index n, focusing on 
resonances based on the stable four-bounce ( "diamond" ) 
PO. The numerical method used is a variant of the "S- 
matrix method" for billiards ]l7| , p^[ . If, as in this case, 
the relevant orbit is stable and we neglect leakage, then it 



is possible to construct 1 19 approximate modes which are 
piecewise gaussian on each segment of the PO. From this 
theory one finds that the effective beam waist in each seg- 
ment will scale as A = £,\/nk r 0l where £ is a constant 
dependent only on the stability matrix eigenvectors of 
that particular segment, and k Q is the quantized eigen- 
value of the mode. In Fig. ^ (a) we plot one represen- 
tative quasi-bound mode at n — 1.56; the corresponding 
far-field angular intensity is plotted in Fig. [j]. 




FIG. 2. (a) Field intensity plot (gray-scale) for a diamond 
resonance of the quadrupole at critical incidence for the points 
at 4> w = 0, 7r, calculated numerically at nk r ~ 90, n = 1.56, 
e = 0.1. Note that there is negligible emission from the up- 
per and lower bounce points at <f) w = ±90° because they are 
above the critical angle (b) Chiral counterpart of this exact 
resonance, simulating a gaussian beam (see text). 

In both figures one sees the rapid oscillations due to 
interference, but in Fig. |l| one can see that the maximum 
of the intensity is displaced from (j) = 90° as expected due 
to Fresnel Filtering. To compare the size of the effect to 
our analytic theory for a planar interface it is convenient 
to eliminate the interference by calculating the "chiral" 
resonance shown in Fig. (b). This is the original reso- 
nance with the negative angular momentum components 
projected out, hence mimicking a uni-directional beam. 
When this is plotted in Fig. |l| (dot-dashed line) it gives 
the smooth envelope of the diamond resonance without 
the oscillations. We can regard this chiral resonance as a 
beam incident at <f> w — on the boundary and compare 
it to our planar model. The angle of incidence of the 
"beam" with respect to the tangent plane is 0i 39°, 
and we have chosen n so that 9i = 9 C ; hence naive ray 
optics predicts tangent emission (4> ma x — 90°). From 
Fig. [I] one sees that the resonance emission is peaked at 
4>max ~ 66°, whereas the planar model gives a similar 
envelope slightly shifted with 4> m ax ~ 62°. In evaluating 
the planar model we use A and zq as calculated from the 



gaussian theory of this diamond resonance, hence we have 
no free parameters. The observed difference is likely due 
primarily to the effect of the curvature of the boundary. 

To evaluate systematically the Fresnel Filtering effect, 
we have calculated the farfield peaks of the set of diamond 
resonances while varying the index of refraction, so that 
the critical angle is scanned through the PO incidence 
angle 0i ~ 39°. In order to remain as close as possible 
to our GM with fixed A we have chosen the resonances 
so that nk r is approximately constant. In Fig. Q the 
exact numerical resonance peak is compared to the cal- 
culated 4>max from Eq. (Q) and to the direction predicted 
by Snell's law. Clearly the deviation from Snell's law, 
A6ff, varies with distance from 9 C ; further studies find 
that this region of significant deviation decreases with 
increasing A as expected. 




FIG. 3. Comparison of peak angular farfield values (f>max 
for varying critical angle 9 C — sin _1 (l/7i). (Diamonds): ex- 
act resonances at nk r ~ 90. (Solid line): GM calcula- 
tion with A ~ 8.82. (Dashed line): Snell's law prediction: 
sin (jyjnax = n sin (9 4 where 9i ~ 39° . A9p F designates the 
deviation from Snell's law at 9 C — 9i. 

In conclusion, we have shown that the transmission 
direction of a narrow beam through a plane dielectric in- 
terface can be quite different from the direction predicted 
by applying Snell's law to the incident beam direction. 
This effect is due to a phenomenon we call Fresnel Filter- 
ing and is of great importance in predicting the emission 
patterns from resonances based on periodic ray orbits in 
micro-cavity lasers. This is true even when the size of the 
resonator ro is much larger than the wavelength and one 
might have expected ray optics to be quite good. Specifi- 
cally the effective beam waist for stable resonances scales 
as A a \/nk r , so from Eq. ([?]) the deviation angle at 
critical incidence 9 C FF tx (rifeo?-,,)" 1 / 4 , and hence may be 
large for nkgro ~ 10 2 — 10 3 as in recent experiments on 
semiconductor ARC lasers • 
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